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Abstract 

Functional equations, in the form of fusion hierarchies, are studied for the trans¬ 
fer matrices of the fused restricted i lattice models of Jimbo, Miwa and Okado. 
Specifically, these equations are solved analytically for the finite-size scaling spec¬ 
tra, central charges and some conformal weights. The results are obtained in terms 
of Rogers dilogarithm and correspond to coset conformal field theories based on 
the affine Lie algebra A^_ x with GKO pair A^}_ 1 © A^f j D A^l 1 . 
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1 Introduction 


1.1 Conformal spectra and lattice models 


In the last decade it has become clear that the critical behavior of two-dimensional 
statistical systems can be described by conformal held theories JT|. |2j. |3|, |J. This is 
possible because statistical systems at critical points possess conformal invariance 0 . In 
particular, the continuum limit of the critical L-state restricted solid-on-solid (RSOS) 
models of Andrews, Baxter and Forrester (ABF) || provide realizations of the unitary 
minimal conformal held theories with central charges 


L{L + 1) 


( 1 . 1 ) 


and conformal weights given by the Kac formula 
[{L + 1 )t - Ls ] 2 - 1 


A-t,S — 


4L(L + 1) 


l<t<L-l, 1 <s<L, s <t . 


( 1 . 2 ) 
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At criticality, the characters of the Virasoro algebra of the conformal field theory appear 
naturally in the modular invariant partition functions. In particular, the conformal 
spectra, including the central charge and the conformal weights, can be obtained from 
the finite-size corrections to the transfer matrix eigenvalues |7], [|. Surprisingly, the same 
characters appear ii in the expressions for the off-critical local state probabilities of 
the ABF models. 

Further conformal held theories are obtained from the fused ABF RSOS models at 
fusion level p. The central charge and conformal weights of these theories are given by 


c = 


At,s — 


6 p 


3 p 

p+ 2 (L + 1 )(L + 1 - p) 

[(L + l)t — (L + 1 - p)s} 2 - p 2 
Ap(L + 1)(L + 1 — p) 


+ 


sp(p~ s 0 ) 

2 p(p + 2) 


where 


0 < so < p and so = ± (t — s) mod 2 p . 


(1.3) 

(1.4) 

(1.5) 


This conformal spectra was conjectured E3 Q by the evaluation of the local state 
probabilities and confirmed by the direct calculation Bi of f ini te-size corrections for 
the fused ABF RSOS models. The fused lattice models are related to coset conformal 
held theories obtained by the Goddard-Kent-Olive (GKO) construction |12|]. The relevant 
GKO pair is 


4 ” 

level L — p — 1 


4 1 * 3 A? 

p L — 1 


( 1 . 6 ) 


The minimal unitary series corresponds to level p — A. For p > 1 the conformal characters 
are identified as branching coefficients. 

Another direction to extend the unitary minimal series is to higher rank theories with 
the GKO pair 

kli-i © d (1.7) 

level l — l 1 l 


These theories reduce to the unitary minimal series for n = 2 and l = L 
Virasoro algebra W" can be constructed 
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1. The extended 
starting from this GKO pair of the affine 
algebra A^ 1 l 1 . It has been confirmed that the integrable lattice models corresponding to 
these higher rank coset conformal field theories (|1.8| ) are the ] lattice models of Jimbo, 
Miwa and Okado JR|. The central charge and conformal weights are 


c = (n — 1) 


n(n + 1) 


(/ + n)(l + n — 1) 


( 1 . 8 ) 


[(/ + n)t — (l + n — l)s] 2 — n(n 2 — 1)/12 
2 (l + n — !)(/ + n) 


(1.9) 
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where t = p + p and s = q + p. Here p is the Weyl vector or the sum of all fundamental 
weights of and p and q are local states given respectively by the dominant integral 
weights P + (n,l — 1) of level l — 1 and P + (n,l) of level /. These notations are briefly 
explained in subsection [2.1| . The conformal spectra (|1.8| ) and o has been obtained by 


the GKO construction of the stress-energy in conformal held theory [16| and the study 
of local state probabilities of the JMO lattice models |T5| . 

Further generalizing ( jl . 7|) and (|1.8| ) by the fusion procedure leads to the GKO pair 


d (1) ffi 4 (1) 4 (1) 

/l n _l ty S±n -1 -> A n-1 


level l — p — 1 


P 


l 


( 1 . 10 ) 


with Virasoro algebra VV"’. Using the GKO construction of the stress-energy tensor, the 
central charge and conformal weights Jl6j] of the W™ models are given by 


c = 


(n 2 — 1 )p n{n 2 — 1 )p 


p + n (/ + n) (/ + n — p) 


( 1 . 11 ) 


^ t,s ~ 


[(l + n)t — (/ — n — p)s] 2 — p 2 n{n 2 — 1)/12 
2 p(l + n — p)(l + n) 


+ h 


v ’ 


( 1 . 12 ) 


where t = p + p, s = q + p and p, q and u are local states respectively in the dominant 
integral weights P + (n,l —p), P + (n,l ) and P + (n,p). By the Feigin-Fuchs construction 
[2j , 5 U is fixed by 


= 


2 pu ■ p — nis 2 
2 p(p + n ) 


(1.13) 


In this paper we calculate the finite-size corrections to the eigenvalue spectra of the 
transfer matrices of the fused JMO ] lattice models. Specifically, we generalize the 
analytic study introduced in H |19[ to obtain the central charges and conformal weights 


of the fused JMO lattice models. The functional equations of the fused transfer ma¬ 
trices for the JMO A^l_ ] lattice model have been given in [17] 


In [20] the thermodynamic 


Bethe ansatz-likc equations (also called y-systems or inversion identity hierarchies) of the 
model have been introduced and by solving these functional equations Kuniba, Nakan- 
ishi and Suzuki have extracted the central charge (|1.11|) from the finite-size corrections 
to the eigenvalue spectra of the transfer matrices. In fact the central charge and con¬ 
formal weights together appear in the finite-size corrections to the eigenvalue spectra of 
the transfer matrices. To obtain the conformal weights we have to consider the excited 
states of the transfer matrices and the calculations are more complicated because of the 
more complicated analyticity for the states. Therefore we generalize the study done in 
[|I8|, [To! f° r the fused JMO A^li lattice models to calculate some conformal weights. In 
the scaling limit the fused JMO models have been shown [17] to yield the same central 


charges ( |1.11| ) as the W” conformal held theories. Here we find the conformal weights 


At,s — 


n(n 2 — 1) [(/ + n)t — (/ + n — p)s ] 2 — p 2 2 pu 


24 


p{l + n — p)(l + n ) 


+ 


nv~ 


2 p(p + n ) 


(1.14) 
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where 1 < s < 


— I 1 

»-iJ’ 


< t < 


n -1 


the largest less than or equal to x. Clear 


and v = 


(s — t) — p , where |_^J denotes 


y, (|1.12|) and (|1.14|) agree if 


(t,s) = 


r = 


s" = 


V = 


(*AP> = 


n(n 2 — 1) 



12 


n[ 

n 2 — 

1) 


12 


n[ 

[n 2 — 

1) 


12 


n[ 

[n 2 — 

1) 


12 


n[ 

[n 2 — 

1) 


is , 


t 2 , 


(1.15) 


12 


/y 


/y 


This conhrms that the underlying solvable statistical mechanics models corresponding 
to the conformal held theories with extended algebra W™ are precisely the fused critical 

[ lattice models of Jimbo, Miwa and Okado. 

The paper is organized as follows. In the next subsection, we describe the fini te-size 
corrections to the eigenvalues of the row transfer matrices of critical lattice models as 
predicted by conformal invariance. In section 2 we define the lattice models of 

Jimbo, Miwa and Okado. In particular, we discuss the fusion rules satisfied by the ad¬ 
jacency matrices of these models and the corresponding functional equations of fused 
transfer matrices. These functional equations can be converted into inversion identity 
hierarchies. These inversion identity hierarchies and their original functional equations 
are the elementary relations needed to determine the finite-size corrections to the eigen¬ 
values of the row transfer matrices. Next we find the asymptotic and bulk behavior of 
these inversion identity hierarchies, which are the important data for obtaining finite-size 
corrections. In section 3, we convert these functional equations into nonlinear integral 
equations which we solve analytically. After using some known |23|] Rogers dilogarithm 
identities, we obtain the finite-size corrections, the central charges and conformal weights. 
Finally, a brief conclusion is given in section 4. 


1.2 Finite-size corrections 

Much work has been done on extracting the conformal spectra of exactly solvable lattice 


models from finite-size corrections. These extensive calculations |24|, |2^, [26], |27|, [28 
BH, 0, pi], |17l |32, I] give very strong and direct evidence to support the predictions of 


conformal and modular invariance. 

According to the predictions of conformal and modular invariance, the partition func¬ 
tion of a two-dimensional lattice model on a finite M x N periodic lattice or torus, can 
be written for large M and N as 


Z m ,n = Tr T m ~ 


Z(q)e~ MNf . (1.16) 

Here T is the row transfer matrix, / is the bulk free energy and Z(q) is the universal 
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finite-size partition function with modular parameter 


q = e 27riT , t = — e ^- hu) . (1.17) 

Suppose that the eigenvalues of the row transfer matrix T of a periodic row of N faces 
are given by 

An = e ~ E " (1.18) 


where E n are the corresponding energy levels. It then follows that 

Z M , N = Y,e~ MEn - (1-19) 

n 


Conformal invariance now 


predicts |33], Q that for large N the energy levels take the 


form 


7TC 

E 0 ~ Nf - — sin (hu) , 

2n 

E n ~ E 0 + — [xn sin (hu) + is n cos {hu)\ 


( 1 . 20 ) 

( 1 . 21 ) 


where E 0 is the ground-state energy, 


x n = A + A + k + k and s n = A — A + k — k 


( 1 . 22 ) 


are respectively the scaling dimensions and spins of the various levels and k, k are integers. 
The numbers c and (A, A) are identified as the central charge and conformal weights of 
the primary operators of the underlying conformal held theory. 


2 The models and their fusion hierarchies 

2.1 Algebraic notations and Boltzmann weights 

A vector a represents a level / dominant integral weight of if 

n— 1 

a = J2 G z + (2.i) 

fi =0 

and E()=o = ^ where Z + is a set of all non-negative integers and A M with /i = 

0,1, • • •, n — 1 are the fundamental weights of , with A„ = Ao- Fix an integer l > 1, 
and denote by P+(n, l) the set of dominant weights. Then an element of P+(n, l ) is called 
a local state. 

An ordered pair of local states (a, b ), with a, b e P+{n, l), is called admissible if 

b — a + jl, fi — 0,1, • • •, n — 1 (2.2) 

where ji are the elementary vectors defined by 

/i = A /t+1 - A ;i ji = 0,1, • • • ,n - 1 . (2.3) 
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Figure 1: The set of local states P + ( 3,4) (n = 3,1 = 4). Each state is associated with a 
corresponding Young diagram and the admissible pairs of the states (a, b ) are represented 
by arrows from a to b. 

For each such pair we place an arrow from a to b, then all local states in set P+(n, l) can be 
represented by an oriented graph. There is a one-to-one correspondence between a state 
in P+(n, l ) and a Young diagram (/i, • • •, /„_i) with l = f 0 > fi > ■ ■ ■ > f n -1 > fn = 0 
(see Figures [L] and |2|) given by 


71—1 

U+ i) a m- 

It is usual to omit the columns of length n in these Young diagrams. 

The lattice models are described by assigning a Boltzmann weight 


W 




to each configuration a, b,c,d 6 P + (n, l ) of four sites surrounding a face. The face weights 
are nonzero only if ( d,a ), (a, b), ( d,c ), (c,b) are all admissible. The nonzero weights are 
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Figure 2: The set of local states P + (4, 2) (n = 4, Z = 2). Each state is associated with a 
corresponding Young diagram and the admissible pairs of the states (a, b ) are represented 
by the arrows from a to b. 


given by 


[35 


W 


a 

a + jl 


a + a 
a + jl + z> 



/t(A + it) 

“MA T’ 

hiXa^ — u ) 
h(Xa^) ’ 


ifO<a = /i = z/<n — 1; 
ifO<a = /i^z/<n — 1; (2.4) 


/i(m) h{\a^ v + A) 
. /i(A) /r(Aa^) 


if 0 < a 


v 7 ^ / j , < n — 1, 


where a is a local state and is given by the inner-product (a + p, fi — u). Here u is 

the spectral parameter, A = p is the sum over all fundamental weights and h[u) 

is the elliptic theta function 


h(u) = sin (it) JJ[(1 — 2 q 2j cos2 u + q 4j )( 1 — q 2j )\- (2.5) 

3 =i 


In this paper we will only consider the critical case when the elliptic nome vanishes q — 0 
and the theta function reduces to the trigonometric function 


h(u) = sin (it). 


( 2 . 6 ) 
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2.2 su(n) fusion rules 

The oriented graphs, such as that shown in Figure [I] and [|. describe all admissible pairs 
(a, b) in P + (n, Z). The corresponding adjacency matrix A with elements 

^ _ f 1 (a, 6) admissible ^ y\ 

a,b ~ 10 otherwise 1 ' ’ 

satisfies certain su(n ) fusion rules which determine the admissible pairs (a, b ) in P+(n, l ) 
for the fused , lattice models as we will now explain. 

For fixed rank n > 2 and level Z, the su(n) fusion rules determine a hierarchy of com¬ 
muting adjacency matrices labeled by the representations of su(n). Set uj\ — (1, 0, • ■ ■, 0) 

and consider the tensor product of two irreducible representations of of su{n ) with 

Young tableaux / = (/j, / 2 , ■ • •, f n -i) and a>i. Suppose the decomposition of the tensor 
product gives s irreducible representations with Young tableaux f k = (/f, f k , • • •, f k _ ,) 
k — 1, 2, ■ • • s with s < n. Then the su{n) fusion rule is expressed as 

A^A = Y j A^ k \ (2.8) 

fc=i 


where = I if / = (0, 0, • • •, 0) and A^ = 0 unless l > fi >■■■ > f n - i > f n — 0. In 
the oriented adjacency graphs s is just the number of outgoing arrows originating from 
the local state with Young tableau / and f 1 , f 2 , ■ ■ ■ f s are the corresponding adjacent 
local states to which the outgoing arrows point (e.g. see Figure |l|). For example, suppose 
that n — 3 and A^ is represented by the Young diagram / = (/j, / 2 , • • •, f n -i)- Then 
the fusion rule (T8) takes the form 




Q 


P 



Q 



t~i ® n 



l : 
^ j 


y n 


9-1 





q +1 


P-1 

n 



9+1 






q 


p +i 





q 



where p — fi — f 2 and q = / 2 . There is a level-rank duality of the adjacency matrices. 
Let / be the Young tableaux obtained from / by transposing the diagram. Then the 
adjacency_matrix A^ of models with rank n and level l is the same as the adjacency 
matrix A^ of models with rank Z + 1 and level n — 1. 


Of particular interest here are weights pA a in P+(n, l) with rectangular Young 
tableaux / = (/ 1 , / 2 , • • •, f n -i) where f s < a = p and f s>a = 0. Let us write the cor¬ 
responding adjacency matrix as Then, from the fusion rule (|2.8|), it follows that 


y\(a,p)^(a,p) _ j^{a,p+l) j^(a-l,p) j^{a+l,p) ^ 9 ) 
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where a = 1,2, 


, n 


1 and p = 1,2,•••,/. Here A^'^ = lifa = 0orp = 0 and the 


closure condition is that A( a,p ' ) = 0 if p < 0 or p > l or a < 0 or a > n. In the next 
subsection [2.3| we will see that a similar relation holds at the level of the transfer matrices 
of the fused models. 


In general, the admissible states of adjacent sites of the fused models are given by 
the su(n) fusion rule. In the appendix we give the explicit solution of the fusion rules for 
n — 4 and 1 = 2. In contrast to fusing the ABF models, the elements of ' for n > 2 
can in general be nonnegative integers greater than one, for example, kb 3,1 ) for n = 3 
and 1 = 5. In such cases we have to distinguish the oriented edges between the adjacency 
states (a, b ) by bond variables a = 1,2,, A 


2.3 Functional equations 


A class of fused A^l j lattice models with adjacency matrices can be constructed by 
the fusion procedure [[I6J. Let us consider the fused models obtained by fusing (a p ) x (b q ) 
elementary blocks corresponding to the representations pA a and qAj,. The fused face 
weights of these models vanish, that is, 


W, 


(M 

(“.p) 


d 

V 

c 

\ 


a 


P 

u = a 

u 

\ a 

li 

b 

) a 

u 


c 

P = 0 
b 


unless A^f ± 0, A^ p) ± 0, A^f ± 0 and A%f ± 0 and 1 < a < A^f, 1 < (3 < A ( °f\ 

1 < p < A and 1 < v < Aj’f\ Suppose that a(a ) and b((3) are allowed state (bond) 
conhgurations of two consecutive rows of a lattice with N columns and periodic boundary 
conditions. The elements of the fused row transfer matrices Tare given by 




N 


u)\ b ,p) = uEK 

i=HvA 


0 b,g) 

( a »P) 


^ &j+i Vj+i bj- i-i 

Pi 


a,- 


V ( ‘j Vi b .i 


\ a*¥. 


U \ ~ 

u 

/ aj 



bj+i 

Pj 


( 2 . 10 ) 


where ajv + i = a 1: b N+ i = bi, Pn+i = r h an d A^ = 0 mod n. Since the fused face weights 
satisfy the Yang-Baxter relation, we obtain a hierarchy of commuting families of transfer 
matrices. Specifically, if (a,p) are held fixed 

[T (6 - 9) (u),T (ft, ’ 9 ')(v)] = 0, (2.11) 


where we have suppressed the subscripts = TP’ q \(u). These transfer matrices 

satisfy a group of functional equations which are expressed in determinant form in [17 


From this determinant form the following useful functional equations can be extracted 

m- 


T {b ’ q \u) T {b ’ q \u - A) = 

T^ q+1 \u) T (6 ’ 9_1) (m - A) + T {b+hq \u ) T^ b ~ 1,q \u - A) 


( 2 . 12 ) 
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These functional equations can be established using the fusion procedure |>b[ and they 
are the Baxterization of the su(n) fusion rule (|2.9|) for the adjacency matrices. 

The functional equations are easily converted into the so called thermodynamic Bethe 
ansatz-like (TBA-like) equations in the sense of |2T| |37|]. Thus inserting 


t^ q) (u) := 


T^ q+1 \u) - A) 


T (b+1 ’ q \u) T (6_1 ’ 9) (w-A) 
into ([2.12|) yields the TBA-like hierarchy 

[I + *(M+i)(„)][! + ^-^(u - A)] 


t M (u)t M (u- A) = 


[I + ( t (ft+h?)( w ))-l][I + ( t (6-l-<z)( M - A))" 1 ] 

The closure condition becomes 

t^°\u) = t {b ’ l) (u) = 0 , b = 1,2, •••,?*- 1 . 

By definition we have the following symmetry. Applying the replacement 

1 


t (M ( 


u 


t u = 




(2.13) 


(2.14) 


(2.15) 


(2.16) 


to the TBA-like hierarchy, we have 

A) = 


i(M i \i(M/ 
t (u)t (u 


[I + i <i,,+1, (u)][I + i ( ‘’ ,1, (ti - A)) 


[I+(t < * + '- ,) ( l .))-'][I+(t ( ' , - 1 - s) ( l <-A))->] 


(2.17) 


with the closure condition 
i(b, o) 


t'"'"'(i.) = t <M (i.) = 0, b = 1,2,.. -,i-l. 


(2.18) 


This symmetry is just level-rank duality 


In terms of the t matrices, the functional equations (|2.12|) can be rewritten as 

T^ p \u) T {a ’ p \u~ A) = T (a ~ lj,) (u-A) T (a+1 ’ p \u)(l + t (a ' p \u)j (2.19) 


It is obvious to see that we need to solve the TBA-like equations ( |2.14|) first to get the 
solutions for the transfer matrices T ta ’ p> (u). 


2.4 Asymptotics of t 

The functional equations ( |2.14|) are easily solved in the braid limit u —> ±ioo. The 
asymptotics t^ q) = t^ b,q \±ioo) satisfy 


+{b,q).(b,q) 


[i + t£^ 1 >][i + tS»- 1 )] 

[i+(«£ +M y i ][i+(i£ _i ' , ' ) )- 1 ] 


( 2 . 20 ) 
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with the closure condition 


+0,o) 

L OC 



b = 1 , 2, • • • , n — 1 . 


( 2 . 21 ) 


Let us write ffi’ 1 ) as 


(ii) = sin [(n + 1)0] 
00 sin[(n - 1)9} 


( 2 . 22 ) 


with 9 to be determined. Using ( |2.2(J|) as a recursion relation for we then derive 

(2.23) 


Jb, q ) = sin[(w + g)g] sin[gd] 
00 sin[(n — b)9\ sin [b9] 


for b — 1 , 2, • • • , n — 1 and q — 1,2, 
quantization 


, l — 1. The closure condition (|2.21|) imposes the 


e = Hhl 


rrij = 1, 2, • • •, 


n — 1 


(2.24) 


where h = n + /. Mathematically, the exponent rrij can go from 1 to h — 1. But t < ^ g ' > 
is on longer the right solutions of the models. This can be checked by 
looking the adjacency matrices or the transfer matrices for small size. 


with rrij > 


We would like to mention here that the solutions given by ( |2.23| ) are not the general 
ones of the equations (|2.20|) . However, the general solutions involve complex eigenvalues, 
which together with the real eigenvalues (|2.23|) give all of the low-lying excited states. 


2.5 Zeros and poles of the transfer matrices 


To find the finite-size correction we need to solve the TBA equations (|2.14|) or (|2.17|) 
in certain analytic domains. Inside these analyticity strips the ground state eigenvalues 
T^ b ' q \u) = T(^p( (u) should not possess any zero except those which come from the 
parameterization of Boltzmann weights. They are of order N and their locations are 
independent of the eigenvalue under consideration. Using self explanatory notation, the 
location of these zeros are as follows: 


uuu m+k-j-m uuand &>« (2.25) 

j=0 1=0 i= 0 k= 0 j=0 1=0 i= 0 

UU U U^ + k-j-W U U U{(?+*-^'-0A} , Q>P and b < a (2.26) 

j= 0 1=0 i= 0 k =0 j =1 1=0i =0 

UU U U U U{( fc +* _ U A } > q<P and b > a (2.27) 

j=0 1=1 i=0 k=0 i=l k=0 j=0 

UU U U{(*+ fc “i“0 A }ij U U{( fc +*“i) A }. Q<P and b < a (2.28) 

j=0 1=1 i=0 k=0 i=0 k=0 j=0 
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These can be seen from the fusion procedure of the models |3S|. To locate the zeros and 
poles of we use the definition ( |2. 13|) and the zeros of T&vju). We distinguish several 
cases as follows: 


(I) 1 < q < p — 1 and b = a 


(II) 


(HI) 


(IV) 


(V) 


(VI) 


zeros: 0 
g 

poles: (J{ZA} 

i=i 

q = p and 1 < b < a — 1 

b 

zeros: LK(i-“W 

3 = 1 

poles: 0 

q = p and b = a 

a— 1 

zeros: U(-M 

J-0 

p 

poles: { ZA } 

i=i 

q > p + 1 and b = a 
zeros: 0 

p 

poles: |J{(/ + g-p)A} 

i=i 

q = p and b > a + 1 

a 

zeros: U((j-a)A} 

3 = 1 

poles: 0 

q 7 ^ p and b ^ a 
zeros: 0 
poles: 0 


(2.29) 


(2.30) 


(2.31) 


(2.32) 


(2.33) 


(2.34) 


This pattern of zeros and poles is divided according to the fusion level (a,p) in the 
vertical direction. We fix the vertical fusion level (a,p) and vary the fusion level (6, q) in 
the horizontal direction. The transfer matrices are free of zeros for q ^ p and free 
of poles for b ^ a. 


2.6 Bulk behavior and the largest eigenvalues 

According to section |2.5| the analyticity strip for t^ a,p ^(u) contains a zero of order N at 
u — 0 and a poles of order N at u = A. All other functions t^ b,q > are analytic and non-zero 
in their analyticity strips — |A < u < |A. For large N the leading bulk behavior to t^ b ’ q ' 
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is given by 


tbtfkiu) = 


constant , 


constant 


tan {\hu) 


q p or b 7^ q , 

, q = p and b = a . 


(2.35) 


The ansatz of bulk behavior matches the zero and pole distribution. The constants are 
fixed by the TBA-like equations (|2.14|) and can be calculated similarly to the asymptotics 
of t( b,q \ Corresponding to the asymptotics solutions (|2.23|) we find that the bulk values 


4uik for 1 < q < p 


1 and 1 < b < n — 1 are given by 


AM _ 

‘'bulk — 


sin[(g + n)a] sin(gcr) 
sin(6cr) sin[(n — b)o\ 


with 


cr = 


p + n 


mb = 1 , 2 , 


P 


n 


and for p + 1 < q < h — n — 1 and l<fe<n — lby 

,(b,q) _ sin[(g -p + n)r ] sin[(q - p)r] 


6 bulk — 


sin(hr) sin[(n — b)r] 


with 


(2.36) 


(2.37) 


(2.38) 


r = 


m j 71 

h — p 


rrij = 1 , 2 , 


l — p 
n — 1 


(2.39) 


Here we see that p — 1,2, ■ ■ ■, h — 2. For the largest eigenvalue, the appropriate choices 
are 9 = n/h, a = n/(p + n) and r = ir/(h — p) in Q2.23|) , (|2.36|) and ( 2.38|) respectively. 
Here for the same reason as the braid limit (|2.24|) we restrict the range of the exponents 




3 Finite-size corrections 

The finite-size corrections for the eigenvalues TAw) can be obtained by solving the equa¬ 
tion 

T m (u)T m (u - A) = T (a-1 ’ p) (w - A)T (a+1 ’ p) (u)(l + f (a ’ p) («)) (3.1) 

Although the JMO A^}_ { models satisfy an inversion relation, there is no crossing sym¬ 
metry for the face weights of these models so we cannot extract the bulk behavior for 
the transfer matrices from the inversion relation alone. Instead the f ini te-size correction 
to T (a w) (u) is given by 


T^p\u) 


rji(a,p) 
1 bulk 


(n)T f 


(a,p) 

finite 


(«)■ 


(3.2) 
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(3.3) 


Inserting (|3.2|) into (|3.1|) and setting the bulk behavior 


as 


Tr 


(ArO, 


bulk 


u)T^(u - A) = T^ a ~ l ' p \u - A )T( a+1 ’ p \u) 


we hnd 


T, 


(“A , 


finite 


u)TtSl(u - A) = 1 + «<“■">(«) 


(3.4) 


We can check that H3.3D is correct for the case of n = 2 or the (l + l)-state ABF models 
in ||. The fini te-size corrections for T^ a,p \u ), therefore, are represented by the inversion 
identity hierarchy t^ a)P \u). In the analytical treatment of (|3.4| ) and (|2.14| ), we will see 
that the only inputs for the finite-size corrections in the scaling limit are the asymptotics 
and bulk behavior. 


3.1 Nonlinear equations for finite-size corrections 


It is useful to introduce functions of a real variable by restricting the eigenvalue functions 
to certain lines in the complex plane, 


il (M l 

{%) 

:= l + a (b ’ 9) (^) , 

a M l 

[x) 

t^( n [x - 



\2h 

b M l 

[x] 

._ T {b ’ q) ( ni a 
finite V2 h 


a — b + p — q 


X 


a — b + p — q + 1 


(3.5) 

(3.6) 

(3.7) 


The functional relation (|3.4|) can then be rewritten in terms of the new functions as 

b M (x - ni/n)b {a ’ p \x + m/n) = U {a ’ p) (a;) . (3.8) 

For the ground state the functions and b (b ’^ are analytic, non-zero in —2n/n < x < 
‘lixjn and possess constant asymptotics for 1Ze x —> ±oo (ANZC). Taking the logarithmic 
derivative of the above equation and introducing Fourier transforms 


B {a ’ p \k) = — / dx [In b (a ' p) (x)]' e~ ikx 
27r J-oo 

/ OO 

dk B {a ’ p \k) e ikx 

-oo 


(3.9) 


with analogous equations for ih“’ p ' and its Fourier transform A^ a,p \ we then find that 


' ' g{n/n)k _|_ g — (7r /n)k 

(3.10) 

Transforming back and defining the kernel k(x) 


k(x ) :=-—-r—r , 

47t cosh(nx/2) 

(3.11) 

we are able to express b ^ a,p ’ in terms of ir a,p \ 


In (, (a ’ p) = k* In U (a,p) + , 

(3.12) 
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where C'- a,p ' > are integration constants. The convolution f * g of two functions / and g is 
defined by 


/ oo roo 

f{x - y)g(y) dy = g(x - y)f(y ) dy . 

-OO J — OO 


(3.13) 


In the case of low-lying excitations we have to take care of the zeros in the analyticity 
strips so that the simple ANZC properties hold. The result ( |3.12| ) is still correct if we 
change the integration path C so that U- a,p >{x) has an ANZC area and Cauchy’s theorem 
can be applied as in ||. The integration constants in (|3.12|) can be evaluated from the 
asymptotics of and b ( ' a ’ p ' > . In this limit (|3.12|) becomes 


lnb£ p) = ±lnil M + . 


(3.14) 


It can be seen that the constants are dependent on the system size N and do not con¬ 
tribute to the A corrections. 


The is from the inversion identity hierarchy and can be solved from (|2.14|) , 

which can be rewritten in terms of ch 6,9 ) as 


a ( b , q ) [x — 'Ki/n ) a ( ' b,q ' > (x + ni/n) ^ (ic)iI (6 ’ 9+1 - ) (a;) 

o( b - 1 >9)(x)a( &+1 -9)(a;) ~~ (x)A (b+1 ' q) (x) ' 

We introduce finite-size correction terms l^ b,q \x ) by writing a^ b,q \x) as 


(3.15) 


a {b ’ q) (x) = 


/(M ( 


x) 


q ^ p or b a 


tanh^ (nx/ 4) l^ a ’ p ^ (x) , q = p and b = a 


(3.16) 


For the ground state all the functions k b ' q \x ) are ANZC in — n < x < n. They satisfy 
the functional equations 


lM ( x _ n i/ n )lM(x + ni/n) (x)il (fe ’ 9+1) (x) 


l(b- l ,q)( x }l( b +l,q)(x) 




(3.17) 


Again applying Fourier transforms to the logarithmic derivative of the equations the 
Fourier transform l/ b,9 )( k) to l^ b,q \x) satishes 


L^ b ~ l ’ q \k) -2cosh(^Ld>’ q \k) + Ld> +1 ' q \k) 

— 1,<?) _J_ J^(b+l,q) _ j^b,q+ 1) _ j^b,q-l) 


(3.18) 


For hxed q we have the closure conditions L^ 0 ’^ = = A^^ = = 0. This set 

of n — 1 linear equations can be rewritten in matrix form as 


(Adj + /i 0 ) • L q (k) = Adj ■ A q - A q+1 - A 9-1 


with 


Jin = — 2 cosh 


/ 7T k\ 

V n ) 


(3.19) 


(3.20) 
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where L q (k ) is [n — 1) x 1 matrix with the elements and the (n — 1) x 1 matrix 
A q (k) has the elements A^ b,q \ The matrix Adj is the same as the adjacency matrix of the 


classical A n _i Dynkin diagram. The equations (|3.19|) are solvable. Similar equations in 
fact appear in || and these can be solved using a similar method. Thus, after integrating 
the equations, we obtain the set of nonlinear integral equations 


la" = U> + K * iU" -k*( /U‘ l+1 + IIP -1 ) + D" , 


(3.21) 


where the entries of the (n — 1) x 1 matrices la q , lU q and lt q are respectively given by 
In a (b ' q) lnil*' 6 ’^ and 


In »“■«>(*) :=(?! q+V °\V a 

y mtanh [nxj 4) , q = p and b = a . 


(3.22) 


The (n — 1) x (n — 1) matrices K and K are the symmetric matrices whose entries in the 
upper right triangle are given by 

1 r°° 

— I dk e‘ n ' J , 

Z7T J — CO \ 


Ki{x) = 


dk e ikxf sh M( n ~ l)vk/n] sinh[j7rfc/n] 


(3.23) 


sinh[7r/c/n] sinh(7rfc) 

K’( X ) = -i r dk ( cothf—j Sinh|( " ~ ). (3.24) 

7r J-oo \ V n J sinh(7T k) J 


The set of equations ( |3.21 ) are also obtained for the excited states after we take care of 
the extra zeros so that the ANZC properties hold in the analyticity strips. 


3.2 Scaling limit 

The f ini te-size corrections can be extracted from the nonlinear integral equations ( |3.21|) 
and ( P.12| ). The system size N enters the nonlinear equations ([3.21 ) through (|3.22|) . The 
function has three asymptotic regimes with transitions in scaling regimes when x is 
of the order of — In N or In N. We suppose that a^ b,q ^ and scale similarly and define: 

e± q \x) := lim e (b,q) (±—(x + In N)] , 

N— >oo \ n J 

a± 9 \x) lim a (b,q ^ (±—(x + lniV)') , (3.25) 

n->oo \ n J 

A± q \x) : = lim mS b,q ^ (±—(x + In N)] = 1 + (a;) , (3.26) 

(V—>oo V n J 


In this scaling limit, (|3.21|) takes the form 

Ha q = £ e g + K * £A q - K * (iA q+1 + iA q ~ l ^j + D q , 

where we suppress the subscripts ±. The entries of the (n-l)xl matrices a 9 , U q and 
t q are respectively given by iA^^ and £e^ b,q ^\ 


(3.27) 


lcS b,q \x) := lna^(a:) , 

£A^ q \x) := In A {b ’ q \x) , 

q p or b a , 

, q = p and b = a , 


&M <*) == {°- 2e - 


(3.28) 

(3.29) 

(3.30) 
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where q = 1, 2, •••,/ — 1 and b — 1, 2, • • •, n — 1 and 

a( b ' q \x) = 0 6 < 1, b > n — 1, g < 1 and q > l — 1 (3.31) 

= 0 6 < 1, b > ii — 1, g < 1 and q > l — 1 (3.32) 


Now let us consider the scaling limit of (/“’^(ic), which gives the finite-size corrections 
to b^ a,p ^(x). Disregarding the integration constants in (|3.12| ), and for fixed x, the finite- 
size corrections to (x) are given by the following expression of order l/N, 

In b ( ' a,p )(x) := lim (k * \niil a ’ p ^){2x/n) 


— lim 

47T N—>oo 


N—> oo 
roc 


dy 


In (y + In N) 


+ 


lniI (o,p ^(—y — In N) 


cosh [re — n(y — In N)/2] cosh [a; + n(y + In N)/ 2] 

/ OO o —^ OO / 1 \ 

d y + — j dy + 0 (-) 


+ 01 


Ntt J- 
cosh x r°° 


Ntt J- 
sinh x f°° 


/_“ e-*(lA^\y) +ML“'" ) M) dy 
/” e-y(tA!^\y) - t.A ( y\y)) dy + o(i) 


m J- c - - — — -- < 333 > 

It can be seen from fl3.27p that for q = p and b = a the integrals in the above equation 

converge because the function £A^ a,p \y) = ln^l + £a^ a,p \y) s j tends to zero faster than 

any exponential. Moreover, as in the case of the critical ABF RSOS models, we have 
M* (y) - £A p _(y) = 0. 

The integrals in (|3.33|) exist and can be evaluated explicitly with the help of the 
dilogarithmic function 

ln(l - y) 


L(x) = — [ dy — -— + ~ lnxln(l — x). 

Jo y 


(3.34) 


To show this, multiplying the derivative of ( |3.16|) with £A^ b,q \ and (|3.16|) itself with 
{£A^ b ' q ^)\ then taking the difference and summing over q and b, and finally integrating 
we derive 

l 1 Ti 1 «00 

EE/ [(£a M )'eA M - £a {b ’ q \£A^ q) y] dx = 

g=l 6=1 

l 1 n 1 r* oo 

EE/ [(^e (M )'M (M - {£e^ q) - D^)(£A^)'] dx , (3.35) 

9=1 6=1 

where the contribution of the kernel cancels due to the symmetry 

k(~x) = k(x) . (3.36) 

Then using the nonlinear integral equations (|3.21|) to simplify the right-hand side and 
integrating the left-hand side of (|3.35|), we can write the f ini te-size corrections as 


l—l n— 1 

EE L 

9=1 6=1 


A (b,q) 


e- y £A {a ’ p \y) dy = 


l—l n—1 

+ IEE D {b,q) £A {b,q) 
9=1 6=1 


(3.37) 
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where the constants D = 0 which can be shown by the asymptotics of the equations 
( p~27p . 

It can be seen from (|3.37|) that the f ini te-size corrections of the transfer matrices in 
the scaling limit depend only on the solutions in the limits x —> ±oo, which are given by 
the asymptotics and bulk behavior. 


3.3 The central charge and the conformal weights 


The following useful dilogarithm identity has been established by Kirillov |23|]. Consider 
the functions 

sin[(g + n)tf>\ sin[g<£>] 


y M (j,n,r) : = 


with 


sin[a<^] sin[(n — a)<p] 
(1 + j)vr 


1 < b < n — 1, 1 < q < r 


V = 


0<j<n + r — 2. 


(3.38) 


(3.39) 


n + r 

It is obvious that they are the asymptotic solutions of the functions equations (|2.14j) with 
r = l or the bulk behavior of the functions equations with r = p and r = l — p. Then the 
following dilogarithmic function identity holds, 

71—1 V -I 

;) 


s(j,n,l ) := L(- ,, . 

4 + y( k ’ m \],n,r)' 

,2 iV „ 


= i(i)( 


(n 2 — 1 )r n(n 2 — 1 )j(j + 2) 


n + r 


n + r 


(3.40) 


6 j 


n 


L2 


n + 1 


+ 6Z_ 


Here the brackets |ycj denotes the greatest integer less than or equal to x. 

The fini te-size corrections (|3.33|) are expressed in terms of the dilogarithm function 

in (ET57D 


pri ch T l -1 n -1 

lnb^(x) 


Nti 


q= 1 6=1 




(3.41) 


The inputs are the asymptotic and bulk solutions obtained in sections [2.4| and [2 . 6[ . Thus 
we have 


In b^ a,p \x) = 


cosh x 
Nti 

cosh x 
Nn 

cosh x 


' n —1 /—I 

EE L 

, 6=1 9=1 

f n -1 l 

EE L 

,6=1 9=1 

^ n —1 


( 1 ^ 

71 — 1 / — I 

-EE L ( 

6=1 9=1 

' 1 )} 

\H( ft ’9)(—oo) / 

,h( 6 ’9)(oo)/ y 

( 1 ^ 

71—1 Z 

-EE L ( 

6=1 9=1 

' 1 ^ 

VH(6.«)(-Oo)/ 

vA( b ’«\oo)J ) 


Nn 


E 


6=1 
n—1 l 


l—p—1 

E ^ 

9=1 


E E ^ 

6 = 19 =p+l 


Ai 6.9) ( — oo) 

1 


+ T 


A(M(-oo) 


71 — 1 l 

EE L 

6=1 9=1 


H(6>p)(— oo) 

1 


H(6-9)(cX)) 


(3.42) 
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With the help of the dilogarithm identity the finite-size corrections can be expressed as 


cosh x 


lnb^“’ p )(x) = ——— (s(m',n,p) + s(m",n,h — n — p) — s(m,n,h — n) ) 


Nn 


7rcoshx /(n 2 — 1 )(h — n — p) n(n 2 — l)m // (m" + 2) 


6N 


h — p 


h — p 


+ 


(n 2 — 1 )p n(n 2 — 1 )m'(m' + 2) 


n + p 


n + p 


(n 2 — 1 )(h — p) n(n 2 — 1 )m(m + 2) 


h 

+ 6 (in' + m" — m) 


h 


n 


n + 1 


+ 6Z . 


(3.43) 


For the excited states in", in', in are greater than 1 and also no longer independent. We 
suppose that 


in' = in — m" + n 


m — m 


p 


(3.44) 


The relation ( |3.44| ) among m", m', m has been shown for n — 2 in [Q and for n > 2 in 
16| . Using this relation and comparing ( p.43| ) with the finite-size correction 


7r ( 


In b p (x ) = —— ( c — 12(A + A) + (k + k)) cosh x + o — , 


6N\ 

(|3.43|) yields the central charge 


N 


(■ n 2 — 1 )p n(n 2 — 1 )p 


c — 

n +p h(h — p) 

and the conformal weights 

n(n 2 — 1 ){[ht — (h — p)s} 2 — p 2 } n(n 2 — 1 )v(nv — 2 p) 
24 ph(h — p) 24 p(p + n ) 

where the exponents s = in and t = m" are integers satisfying 


1 < s < 


l 


n — 1 


1 < t < 


l — p 
n — 1 


, p= l,2,---,h-2 


and 


v := (s — t) — 


s — t 

p 


p 


(3.45) 


(3.46) 


(3.47) 


(3.48) 


(3.49) 


The integers k and k are such that 


k + k = n(n 2 — 1) 
= 0 mod 6 


s — t 

p 


n 


s — t 

p 


+ 2 + 6n 


s — t 

p 


n + 1 


+ 6Z 


(3.50) 


do not contribute to the central charge and conformal weights and hence can be discarded. 
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4 Discussion 


We have presented the results for the finite-size corrections of transfer matrices for fused 
A^li models of Jimbo, Miwa and Okado. These calculations generalize the evaluation of 
the finite-size corrections of transfer matrices for the fused ABF models given in || and 
the central charge calculation of the models given in ||19|| . In this paper we have 

studied the fused transfer matrices with rectangular Young diagrams 


[a,p\ 



These transfer matrices satisfy the functional equations ( 2.12| ). The finite-size corrections 
of the transfer matrices with the fusion of Young diagram [a,p\ in vertical direction are 
expressed by (|3.43|) . Using Young diagrams the result ( [3.43|) can be represented by 


^2 ([b, l + n - p] + [b, P + n] - [ b , l + n]\ 

6=1 k ' 


(4.51) 


Suppose that the solution ( |3.38| ) corresponds to the Young diagram [b, q] . The expression 
(|4.51|) gives the Rogers dilogarithm explanation for the central charges ( |1.11|) and the 
conformal weights (|1.14|) , which are missed in [^3|. The expression ( |4.51|) also shows 
the coset structure corresponding to the central charge ( |1.11| ) and the conformal weights 
(|1.12|) with the case of ( 1.15 ). For other cases the similar coset has been conjectured [f[(] 
and need to be confirmed, which have not been done in this paper. 

The conformal weights obtained in this paper are a subset of all possible conformal 
weights (|1.12|) . But our calculation has shown that in order to have all conformal weights 
we have to analyze the complex eigenvalues of the inversion identity hierarchies of the 
models. The difficulty part is to find the related dilogarithm identity involving the 
complex eigenvalues. These interesting questions are left for further investigation. 

The method by solving the functional equations of fusion hierarchies to fold the finite 
size corrections of the eigenvalues of the fused transfer matrices firstly was developed in 
for study of the fused ABF restricted SOS models and more recently for the central 

In fact, the functional equations of the fused 


charges of the fused ] models ||T9 


transfer matrices contain enough information to calculate the eigenvalues of the transfer 
matrices. Hence the eigenvalues and the relevant Bethe ansatz equations can be extracted 
from the functional equations and the finite size corrections can be calculated by solving 
the Bethe ansatz equations 00- These two methods technically are different, solving 
Bethe ansatz one employs the string hypothesis and the other one relies on the ANZC 
property. But it seems that by solving the functional equations it is more powerful to find 
the conformal weights. Our present paper has shown this partially for the fused ] 
models. There are a large group of fused transfer matrices for the higher rank n > 2. For 
example, the fused transfer matrices with other Young diagrams. These have not been 
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evaluated here. But the finite-size corrections of these transfer matrices should fit the 


picture (|4.51|) . However, we need other functional equations to solve the fused transfer 
matrices with other Young diagrams. 

There are many solvable IRF models, e.g. models in |I(J, Aff models 

in |Hj and dilute Ai models in HU- Also we have the functional equations for the fused 
transfer matrices of the models [|T^, [HI]. ^3|, [P], [To||. It should be possible to solve these 
functional equations by similar methods. 
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Appendix A: The adjacency matrices of the fused JMO , 
models for n = 4 and l = 2. 

The set of dominant integrable weights P + (4, 2) is given by [(0,0,0),(1,0,0),(0,1,0), 
(0,0,1),(2,0,0),(1,1,0),(1,0,1),(0,2,0),(0,1,1),(0,0,2)] according to the Young diagram no¬ 
tation (see Figure |2j). The adjacency matrices of the model with rank n = 4 and level 
l = 2 is labeled by the order of the elements in P+(4, 2) for the rows and the columns of 
matrix. 
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